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We study decoherence induced by a dynamic environment undergoing a quantum phase transition. 
Environment's susceptibility to perturbations - and, consequently, efficiency of decoherence - is 
amplified near a critical point. Over and above this near-critical susceptibility increase, we show 
that decoherence is dramatically enhanced by non-equilibrium critical dynamics of the environment. 
We derive a simple expression relating decoherence to the universal critical exponents exhibiting 
deep connections with the theory of topological defect creation in non-equilibrium phase transitions. 
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I. INTRODUCTION 

Decoherence of a qubit coupled to an environment is 
studied with the central spin models. In addition to be- 
ing a workhorse for understanding of quantum to classi- 
cal transition p]-[3|, central spin models describe loss of 
qubit's coherence in NV centers in diamond Q, quan- 
tum dots in semiconductors [5| , NMR experiments p, [7j , 
etc. The environment, usually regarded as a destroyer of 
quantum coherence and entanglement, has been recently 
proposed to enhance sensitivity of a magnetometer based 
on the central spin system Therefore, understanding 
how the environment influences qubit's state is central to 
quantum computing and metrology. 

So far, studies of decoherence focused mainly on sys- 
tems prepared in Schrodinger cat-like superpositions in- 
stantaneously coupled to an otherwise static environment 
described by time-independent Hamiltonian. We con- 
sider instead a spin- 1/2 system coupled to a driven envi- 
ronment. The environment undergoes a quantum phase 
transition (QPT) as it is driven through the quantum 
critical point by the change of the external field (e.g., 
a magnetic field in the Ising chain depicted in Fig. [1]). 
Decoherence occurs because the environment "monitors" 
the system and "finds out" its pointer states (i.e., effec- 
tively classical states preserved in spite of decoherence 
[1-3]). This means that each pointer state remains un- 
touched, but makes a distinct imprint on the quantum 
state of the environment. How close to orthogonal are 
these record states determines how well the environment 
knows about the system. Their overlap is the decoher- 
ence factor d(t) multiplying off-diagonal terms of the re- 
duced density matrix of the qubit. 

Effectiveness of decoherence depends on how easy it 
is to perturb - "write on" - the environment: Its "stiff- 
ness" is determined by its Hamiltonian as well as by its 
state. For the Ising environment, Hamiltonian-induced 
stiffness is due to either spin-spin interactions or the ex- 
ternal field (Fig. [1]). When the two balance, the envi- 
ronment is at the critical point, where it is most sus- 
ceptible to perturbations. Decoherence is expected to be 
strongest there. This was already seen in static critical 
environments 0, [To| . In our system the environment is 
not static: It is driven across the phase transition, so 



its instantaneous state - and thus decoherence - is de- 
termined by both its Hamiltonian and its history. We 
show that the latter profoundly affects qubit's decoher- 
ence rate. 

The aim of this paper is to understand how a driven 
critical environment affects decoherence of a qubit. For 
clarity of our discussion, which can be extended to more 
complicated systems, we consider an exactly solvable 
paradigmatic model of a QPT: the Ising chain coupled 
to a central spin- 1/2 (Fig. [1]). We show analytically that 
decoherence is enhanced by the nonequilibrium transi- 
tion and exhibit interesting dynamics characterized by 
the transition rate, coupling strengths and ground state 
fidelity. We also show that the time-dependent deco- 
herence factor encodes both universal properties of the 
system associated with quantum criticality (critical ex- 
ponents) and rich system-dependent (non-universal) fea- 
tures. Our work combines theory of decoherence with 
rapidly growing fields of dynamics of QPTs [TMl and 
quantum information [l7l-[2Cij. In the latter context, our 
study can be regarded as a dynamic generalization of the 
so-called fidelity approach to QPTs [U [H| . We would 
also like to mention that in the context of dynamics of 
QPTs decoherence has been studied in |2ll-[23|. 



II. BASICS OF QUENCH DYNAMICS 

To study non-equilibrium environment - qubit evolu- 
tion we characterize dynamics of an environment driven 
across a critical point. This is done using the quantum 
version [IM1 of the Kibble-Zurek (KZ) theory HHH. 
The environment is quenched by decreasing the external 
field 



git) 



t/ r Q, 



(1) 



where g c stands for strength of the external field at the 
critical point, the quench time tq is inversely propor- 
tional to the quench speed, and time t goes from — oo 
to +oo. A QPT is characterized by vanishing excita- 
tion gap Ao|g — g c \ zv and divergent correlation length 
£o/|# — 9c\ v i where z and v are the critical exponents, 
while Ao and £o are constants [26j]. There are two rel- 
evant time scales during a quench: the system reaction 
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time given by h/A^g — g c \ zv and the time scale of change 
of the Hamiltonian given by (g(t) - g c )/-^(g(t) - g c )- 
Away from the critical point the reaction time is short 
compared to quench time scale and so evolution is adia- 
batic. Near the critical point opposite happens and evo- 
lution is approximately diabatic. The border between the 
two regimes lies at a distance g from the critical point. 
Equating the two timescales, we get 



-l/(l+zv) 



(2) 



KZ theory also predicts that there is a non-equilibrium 
length scale £_imprinted onto the post-transition state of 
the system [25j. It is given by the coherence length at 
the border between adiabatic and impulse regimes where 
the state of the many-body system freezes out: 



(3) 



This scale yields e.g. a typical distance between topolog- 
ical defects created during a symmetry-breaking QPT. It 
also suggests that non-equilibrium excitations resulting 
from the quench should encode a characteristic momen- 
tum scale 



C 1 ~t q u/{1+zv) . 



III. THE MODEL 
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We consider decoherence of the qubit coupled to a 
driven Ising chain undergoing a QPT (Fig [I]). The Ising 
Hamiltonian is 



N 



while its coupling to the qubit reads 



N 



Above the environment size N ^> 1, the system- 
environment coupling 5 <C 1, and a z s is the Pauli matrix 
of the qubit. Ising chain in a transverse field has been ex- 
perimentally studied in [27], [HI . The coupling Hss has 
been implemented in The total qubit-environment 
Hamiltonian is given by 

H — Hg + Hss • 

The driving is a result of a slow - tq ^> 1 - ramp down 
of a magnetic field 



g(t) = 1 - t/T Q , 



(5) 



which is equivalent to (pp) taken with g c = 1. 

Evolution starts at t = — oo with the qubit in a pure state 

superposition of c+ |t) + c_ |^) and the environment in 




FIG. 1: (color online). The central spin- 1/2 (qubit) interacts 
with the Ising spin environment whose quantum phase tran- 
sition is driven by an external magnetic field. For large initial 
and final fields - right and left panel, respectively - the en- 
vironment is in the paramagnetic phase: its spins align with 
the field. For small fields (\g(t)\ < 1, center panel) Ising chain 
enters a ferromagnetic phase, where its spins try to align with 
x or —x breaking the symmetry of the Ising Hamiltonian. 



its instantaneous ground state \GS(t = — oo)). The com- 
posite wave function of the system is given initially 

by 

\^{t = -oo)) = (c+ It) + c_ |D) ® \GS(t = -oo)). 
A straightforward calculation then shows that 

|V>(*)> = fe- i ^ dt " {9{t)) \^{t = -oo)) 

= c+| t)®Te- i/! ~ (»(*)+«) = -oo)) 

+ c_| l)®fe- iI -~ dt " £i9it) - s ' ) \GS(t = -oo)) 
= c+\\)®\<p + (t)) + c-\l)®\<p-(t)), 

where T is the time-ordering operator, and evolution of 
the environmental states coupled to up-down qubit states 
is given by 

i^\cp ± (t)) = H £ ( 9 (t)±5)\cp ± (t)). (6) 

Therefore, evolution of the system depends on the dy- 
namics of two Ising branches evolving in an effective mag- 
netic field given by g(t) ± 6. 

The reduced density matrix of the qubit in the 
{I t)? I i)} basis reads 



PS (t) = Tr £ \m)(m\ 



|c + | 2 c + c*_d*(t) 
c* + c-d(t) |c_| 2 



where d(t) = {tp+(t)\tp-(t)) is the decoherence factor. We 
study its squared modulus, 

D = \d(t)\ 2 = \(v + (t)\^(t))\ 2 , 

also known as the decoherence factor; we follow this 
nomenclature below. When D = 1, the qubit is in a 
pure state, but when D = it is completely decohered. 
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IV. EXACT SOLUTION 



DYNAMICS OF DECOHERENCE 



Dynamics of the decoherence factor can be obtained 
analytically by solving ([6]). This can be done after 
mapping the spins onto non-interacting fermions via the 
Jordan- Wigner transformation [l4|- We introduce after 
Dziarmaga the Bogolubov modes and through 



\<P± 



(*)> = n io*,o_ fc > - v ± w m,i- fc >] , (7) 



fc>0 



where \mk : m-k) describes the state with m = 0, 1 pairs 
of quasiparticles with momentum k = (2s + l)7r/iV, s = 
0, 1, A/"/ 2 — 1. We consider even N. In this formalism 
one easily finds that 



D(t) = Y[F k (t), F k (t) 

k>0 



«r(*)«fc(*)+«r(*)«fcW 



(8) 

We will refer to Fk to describe dynamics of decoherence 
from the momentum space angle. We call it a decoher- 
ence factor in momentum space. Using (j8]) we find that 
for thermodynamically large environments 



(9) 



The above integral representation shall be accurate for 
time instants t for which In Fk (t) is non-singular Vfc [29j . 
All we need now are the Bogolubov modes at various 
stages of evolution. The modes evolve according to Q 



d 
dt 



dt 



—2(g(t) ± S — cos k)vj: + 2 sin fcfxjjr, 
2(g(t) ±5 — cos k)u^ + 2 sin fc^. 



We apply the transformation t± = 
4tq sin k (— g(t) =F 5 + cos fc). Defining Tq = 4tq sin 2 fc 
we end up with the Landau- Zener system, 



. d 

l dK 




(10) 



that can be solved with the Weber functions (see e.g. 
[30|). The initial conditions are v^(t± = — oo) = and 
uf(t± = — oo) = 1. To express the solution in a compact 

form we introduce: z± = t'± ex.p(iir / 4) / ^Jr^ and n = 

—itq/A. The exact solution reads 



■ exp (-ttt^/16) P_ n _i(^±) (11) 



<® - 

uf(t) = exp (-tttq/16) exp(i7r/4) x 

[(1 + n)V_ n _ 2 (iz±) + iz±D- n - 1 (iz±)] (12) 

where V m (iz) is the Weber function [3l|. It is simplified 
in the Appendix. 



The quench starts with the environment in the param- 
agnetic phase (g(t) > 1), moves the environment across 
the ferromagnetic phase ( — 1 < g(t) < 1), and finally 
brings it to the "other" paramagnetic phase (g(t) < —1). 
The critical points are located at g c = ±1. Substitut- 
ing z,v = 1 into (|2]) and ((3]) results in g ~ and 

£ ~ V^Q 0- ^ ne non-equilibrium length scale £ pro- 
vides the typical spacing of kinks (spin flips) excited dur- 
ing the transition [l2|, [3 • Dynamics of the decoherence 
factor is depicted in Figs. [2JU 

As the field polarizing spins in the chain is ramped 
down ([5]), environment's sensitivity to external influence 
increases. This results in enhanced decoherence visible as 
a gradual decrease of the decoherence factor D away from 
the critical point (Fig. [2]). Initially evolution proceeds 
adiabatically and D is identical to ground state fidelity 
of the environment [13, EH, [20| . Namely, to the squared 
overlap between the ground states of Hg (g±5): | (GS(g + 
5)\GS(g — S)}\ 2 . As is derived in the Appendix, away 
from the critical points, \g(t) — g c \ ^> 5, we have for the 
k modes evolving adiabatically 



F k (t) = 1 



5 2 sin 2 k 



(l-2g(t) cosk + g(t) 



2\2 



+ 0(5% (13) 



which compares well to numerics (Fig. EK; see also Figs. 
[St and [5^). Using (j9]) we obtain in the paramagnetic 
phase that 



D(t) « exp 



NS 2 



±g(t) 2 (g(t) 2 - 1) 



(14) 



Just above the first critical point, 1 < g < 1 + ^, the 
adiabatic approximation breaks down and decoherence 
speeds up; decoherence factor decreases substantially. 
This happens because (i) environment's sensitivity is en- 
hanced by quantum criticality amplifying the perturba- 
tion resulting from the coupling to the qubit system; (ii) 
the environment is no longer in the ground state: quench 
is taking it into a superposition of energy eigenstates. 

As the environment is driven past the first critical point 
at g c = 1, there are either partial revivals (Figs. [2] and 
|3} or a monotonic decay (Fig. [4]) of coherence between 
the critical points. They result from non-equilibrium dy- 
namics around the critical point leading to excitation of 
small k modes (the gap in the excitation spectrum closes 
for k = at g c = 1). Expanding (JTTJ and (fl2|) - see 
the Appendix - we find a remarkably simple and accu- 
rate solution for the problem. For low k modes, i.e., 
k - k = 0(l)/y/TQ Q, we get 



F k (t) 



1 



-2"7ttq/c 



: ) sin 2 (15) 



where (j)(t) ~ 4(1 — g(t))rQS = At5. As shown in Fig. 
O) - see also Fig. [5]i - there is a very good agreement 
between ([T5]) and numerics. The larger k modes, those 
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paramagnetic 



paramagnetic 



g =-i 
magnetic field g(t) 



FIG. 2: (color online) Decoherence during a quench: a quasi- 
periodic regime between the critical points. Black solid line 
shows the result obtained from numerical integration of (|10j) . 
Our analytical approximations are superimposed on it as 
a blue dashed line: {16} for g(t) £ (-1,1) and {19]) for 
g(t) < —1. The green dashed-dotted line is the adiabatic 
result equivalent to ground state fidelity (tq — > oo limit): 
(|A1[) combined with (8]) . Numerical solution almost perfectly 
overlaps with analytical results away from the critical points: 
\g(t) — g c \ > g,S. The main plot shows that for moderate en- 
vironment sizes - N — 1000 here - almost perfect revivals of 
coherence take place between the critical points, but coher- 
ence is virtually lost after driving the environment through 
the second critical point. The inset shows the same, but for 
N — 30, 000. It illustrates that for very large environments 
there is little coherence left after passing the first critical point 
and the qubit is completely decohered after crossing the sec- 
ond critical point. Note that ground state fidelity provides 
the envelope for the revivals of coherence between the critical 
points. We used here S — 10 -2 and tq — 250. 



for which k ^> k, evolve adiabatically through the critical 
point. Their decoherence factor Fk reads (fl3|) . which is 
illustrated in Fig. [5b, where again numerics matches well 
our theory. 

One striking feature of ([T5]) is that there is a charac- 
teristic momentum scale imprinted by the quench. It is 
predicted by the KZ theory (j4]). Indeed, (fT5j) depends on 
momentum through the combination of k/k. Moreover, 
the modes centered around 



ln2 1 



contribute most to non-equilibrium decoherence result- 
ing from crossing the first critical point: ([T5]) has a single 
minimum at k m . Interestingly, the position of this mini- 
mum is fixed in the post-transition state, dk m /dt = for 
g(t) < 1. It means that the characteristic length scale 
imprinted by the non-equilibrium critical dynamics re- 
mains "frozen" after the critical point was passed. This 
is in perfect agreement with the adiabatic-impulse pic- 
ture of dynamics proposed by the quantum KZ theory 




g^" 1 g u 

magnetic field g(t) 

FIG. 3: (color online) Same as in Fig. [21 except we consider 
here N-th root of the decoherence factor D illustrating details 
of dynamics of decoherence in the large N limit. Numerical 
result (black line) is obtained by calculating D 1 ^ for N — 
2000 environmental spins. iV-th root of ground state fidelity 
(green dashed-dotted line) is calculated for the same N. 



mm. 

Combining <jT3j) , (jT5]> . and (j9| we get 



D(t) « exp 



NS 2 



4(1 -g{t) 2 ) 



exp 



N fW)] 
2tt 



, (16) 



where f[<j>] = --j= f ™ ds - 4(e~ s2 - e~ 2s * '') sin 2 0] 

and (j) ~ 4tS (the tq — >• oo limit was assumed to simplify 
the upper integration limit in f [(/)]) [29]. The first factor 
above is again ground state fidelity, but this time in the 
ferromagnetic phase: notice the lack of g~ 2 factor as com- 
pared to the paramagnetic phase result (fH|) . The second 
one comes from non-equilibrium environment excitation 
happening around the first critical point. Indeed, since 
we already know that the quench enforces Fk = F(k/k) 
for low momentum modes, a simple change of the integra- 
tion variable to k/k in (|9]) gives the Tq V ^ 1+zv ^ = Tq 1 ^ 2 
factor in ([T6]) ; see also Sec. EH Thus, we show analyt- 
ically that for the Ising chain non-equilibrium contribu- 
tion to the decoherence rate encodes universal informa- 
tion about the criticality of the environment via the z and 
v exponents. This information is recorded in the decoher- 
ing state of the simplest quantum system, the spin- 1/2 
qubit, and should be experimentally accessible. 

The non-universal contribution to non-equilibrium de- 
coherence, f[(j)(t)\ from ([T6]h leads to rich dynamics 
of decoherence. We observe partial revivals of co- 
herence between the critical points taking place when 
the environment-qubit coupling is strong enough: 5 > 
7t/16tq. They happen in the magnetic field (g) domain 
with the period of tt/Atq5. The envelope of these re- 
vivals is given accurately by ground state fidelity: at the 
peaks of the decoherence factor non-adiabatic Fk's ([T5]) 
equal unity and so only adiabatic modes contribute to 
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-10 -5 
magnetic field g(t) 

FIG. 4: (color online) Decoherence during a quench: a gaus- 
sian regime between the critical points. Solid black line shows 
the result coming from numerical integration of (|1Q|) . The 
dashed blue line shows our analytical approximations: the 
gaussian result f)17|) between the critical points [g(t) G ( — 1,1)] 
and (fT9|) past the second critical point [g(t) < — 1]. The in- 
set shows N-th root of the data from the main plot: it il- 
lustrates details of the revivals of coherence past the second 
critical point. The parameters for this plot are: N = 1000, 
5 — 4 x 10 -4 , and tq — 250. Except for 6, they are the same 
as in Fig. [2j The switch from coherence revivals (Fig. [2]) to 
monotonic decay of coherence between the critical points is a 
result of lowering the qubit - environment coupling S below 
the 7t/16tq threshold. 



the decoherence factor there. As ground state fidelity is 
smaller than unity, the revivals of coherence cannot be 
perfect (Fig. [2]). The presence of the adiabatic "fidelity" 
envelope, however, may allow for measurement of ground 
state fidelity at the peaks of the decoherence factor: an 
interesting observation as we have here an out of equilib- 
rium system and ground state fidelity is an equilibrium 
quantity. 

In the opposite limit of weak coupling between the en- 
vironment and the qubit, S < 7t/16tq, decoherence factor 
decreases monotonically till the second critical point at 
g c = — 1. In particular, when 5 <C tt/Wtq expansion of 
([T6]) gives gaussian decay of D (Fig. 2]): 



D(t) « exp 



-N86H 2 (V2-1) 



exp 



-N5 2 



4(1 -g(m 



(17) 

Even richer dynamics of decoherence emerges when the 
magnetic field reverses polarity and drives the system 
past the second critical point at g c = — 1 (Figs. [2JH]). 
Indeed, at the second critical point the modes k ~ tt — k 
become excited because the gap closes there at k = tt. 
This means that the decoherence factor Fk will be given 
by (HSJ for k ~ fc, CEj) for k < k < tt - fc, and 



F k (t) = 1-4 (e-^TQik-nf _ e -^ TQ {k-^ gin 2 ^ 

(18) 

for k ~ 7r — fc, where rj(t) ~ 4(g(£) + 1)tqS (see the Ap- 
pendix). This is illustrated in Figs. [SJi-f, where good 



agreement between theory and numerics is confirmed. 
Similarily as above, the momentum scale from the KZ 
theory leaves a long lasting imprint here: the modes 
most contributing to non-equilibrium decoherence due to 
crossing of the second critical point are centered around 
k = 7r — k m for g{t) < —1: (fT8]) has a minimum at such 
k. 

Marrying CE}, JED and {TED with ® we obtain [29| 



D{t) 



exp 
exp 



NS 2 



2tt 



f[v(t)] \ 

T Q J 



(19) 



Comparing (fl9|) to (fl6j) we see one additional factor, 
exp(— Nf[rj(t)]/27T^fQ) 1 due to non-equilibrium transi- 
tion across the second critical point. We would like to 
stress again that dynamics across the two critical points 
happens "independently" because the gap of the environ- 
ment closes at g c = 1 (g c = — 1) for ko = (fco = tt): the 
first (second) transition excites small (large) momentum 
modes. Excitations created around the two critical points 
impose the same revival period for D, but their phases 
are shifted. It leads to destructive interference damping 
the amplitude of the partial revivals dramatically. We 
also note again that (jl3til9j) are derived away from the 
critical points, i.e., for \g(t) — g c \ ^> 5, g. 

Finally, we compare the quench- amplified decoherence 
factor to the one resulting from the adiabatic quench 
(tq = oo ). We see in Figs. [2] and [3] that decoherence 
can be much stronger in the driven environment: com- 
pare the green dashed-dotted line to the blue dashed line. 
This is related to the excitations (kinks for the Ising chain 
[l2l [lij ) appearing in a non-adiabatic quench. Thus, 
susceptibility of many-body environments to perturba- 
tions reflects not just their instantaneous Hamiltonians or 
eigenstates, but also their actual (non-equilibrium) state 
encoding the quench history. 



VI. GENERAL SCALING RESULT 

A variety of systems have the decoherence factor D 
given in the product form J5J). This factorization ap- 
pears when a many-body environment is described by a 
product wave function in momentum space such as (0. 
This is the case for some quantum magnets (XY, Ising, 
Kitaev, etcO [TBI, conventional superconductors (either 
electronic [32] or those emulated by cold atoms [33]), and 
liquid 3 He [34|. As for the Ising case, excitation of these 
environments during a dynamical phase transition shall 
imprint a non-equilibrium length scale £ and momentum 
scale fc, (|3j) and dU), respectively. Note that we again as- 
sume that the external field driving the QPT is quenched 
according to (pQ). It is thus natural to assume that for 
those systems as well overlap Fk between dynamically ex- 
cited states of the environment depends on momentum 
through (k — ko)/k combination: Fk = F((k — ko)/k). 
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Before reaching the first critical point: 



1-10' 





Between the critical points: 
1 




1-10 

0.1 0.1 71/2 
After pasing the second critical point: 

1 



momentum k 




1-10' 

0.12 0.12 71/2 3.04 3.04 



FIG. 5: (color online) Decoherence factor Fj~ during the 
quench. Solid black lines show numerics based on (|10p . The 
blue dashed lines illustrate our analytical approximations. 
Before reaching the critical point at g c = 1 modes evolve 
adiabatically: panel (a) made for g(t) = 2. The adiabatic 
result (|13|) fits numerics for all momenta. Between the crit- 
ical points - panels (b) and (c) made for g(t) — - low k 
modes are excited due to driving across g c = 1, while the 
large k modes evolve adiabatically. The former is illustrated 
in panel (b) where JT5) fits numerics, while the latter is de- 
picted in panel (c) where (|13|) matches numerics. Panels (d-f) 
prepared for g(t) = —2 show what happens after passing the 
second critical point at g c = — 1: low k modes are still ex- 
cited due to driving across g c = 1, intermediate modes evolve 
adiabatically, while large k modes are excited by crossing the 
g c — — 1 critical point. This is theoretically described by (|15|) . 
([T3]) and JT5J) depicted in panels (d), (e) and (f), respectively. 
Note that despite the tiny range of variation of the adiabat- 
ically evolving modes - panels (a), (c) and (e) - the striking 
envelope of the quasi-periodic revivals of coherence between 
the critical points comes from these modes. Solid black lines 
(numerics) have been obtained for N = 2000, tq = 250 and 
6 — 10 -2 . The same tq and 8 have been used to get the 
dashed lines (analytical approximations). 



where ko is the momentum for which the gap closes. 
When that happens, exchanging the product ([5]) into an 
integral similarly as in (J9j), one can easily derive that 
contribution to the decoherence factor coming from low 
energy modes excited by the quench is given by 



D(t) « exp (-ATT(t)/r^ /(1+zly) ) . 



(20) 



Above 



rv I {l-\-zv) 



captures the universality-class depen- 



dence of decoherence caused by non-equilibrium crossing 
of the quantum critical point, while Y(t) stands for the 
non-universal (system-dependent) response to driving. 
Moreover, r is the environment dimensionality and TV 
is the environment size (number of spins/particles/etc.) . 
Note also that the contribution of the momentum modes 
that were not excited by the quench adds another factor 
besides D to the full decoherence factor D in complete 



analogy to what we have seen for the Ising environment in 
Eqs. ([T6]) and JTUJ). We expect that this "adiabatic" fac- 
tor is well approximated by instantaneous ground state 
fidelity as was shown to be the case in the Ising model. 

Finally, we mention that ([20]) reveals also that non- 
equilibrium decoherence is expected to be exponential in 
the number of topological defects created by the quench: 

N/t ~N/t™ /{1+zv \ 

Note that £ provides a typical distance between topolog- 
ical defects created by the quench. 



VII. CONCLUSIONS 

We have demonstrated rich manifestations of deco- 
herence induced by a dynamical phase transition in the 
environment. In particular, we have shown that near 
the critical point decoherence is dramatically enhanced 
not just by the increased susceptibility, but also by 
the excitations that depend on the past history of the 
system. Our results provide a strong indication of a 
universality-class dependent behavior that bears remark- 
able similarity to the dynamics of defect formation in 
non-equilibrium phase transitions 0, HH • 

This work should stimulate future theoretical studies. 
In particular, it would be interesting to study both uni- 
versal and system-dependent contributions to decoher- 
ence dynamics in other out-of-equilibrium environments, 
especially those described by a different universality class 
than the Ising model. It would be also interesting to find 
out if a similar description can be worked out for sys- 
tems whose wave function is not given by the product 
over momentum modes. Generalization of our results to 
finite temperatures can provide yet another extension of 
our work. 

We expect that our findings will be experimentally ac- 
cessible in future cold atom/ion quantum simulators of 
condensed matter systems [35] . Moreover, it should be 
also possible to extend them to optomechanical systems 
where a cold atom cloud undergoing a QPT (the environ- 
ment) is coupled to a mirror (a quantum system) [36j . Fi- 
nally, generalizing the idea from ^ , the non-equilibrium 
enhancement of the decoherence rate of the qubit might 
find its applications in magnetometry. 
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Appendix : Simplification of the exact solution 

Exact solutions (fTT|) and (fl2|) for the Bogolubov modes 
can be simplified for (i) modes evolving adiabatically; (ii) 
small momentum modes excited by crossing the critical 
point at g c = 1; and (iii) for large momentum modes ex- 
cited due to non-adiabatic crossing of the critical point at 
g c = —1. The resulting expressions can be used to sim- 
plify decoherence factor in momentum space: Fk. Their 
knowledge allows for accurate determination of the deco- 
herence factor D away from the critical points. We will 
briefly outline below how compact expressions for Fk can 
be derived. 

Adiabatic evolution of modes: Modes evolving adi- 
abatically are described by well-known ground state so- 
lutions taken at instantaneous value of the magnetic field 
g(t) [li,[26|. This implies that 

v± = cos^/2), uf = sin(^/2), 



where ± e [0,7r], cos(6> ± ) = e±/yl + e±, and e± = 

T $ + cos k) I sin k. The decoherence factor in mo- 
mentum space takes the compact form 



Fk = cos 



o+ + o- 



(Al) 



which can be simplified by the Taylor expansion 
5 2 sin 2 k 



F k = l 



(1 - 2g cos/c + g 2 ) 2 



0(5*). 



Small momentum expansion for the modes ex- 
cited after crossing the critical point at g c = 1: We 

assume here that g(t) < 1 and consider only k <C 7r/4 
relevant for slow transitions (tq ^> 1). Our expansion 
of (JTTJ) and (fT2|) is based on large \iz±\ expansion of the 
Weber functions. We define z± = a± exp(z7r/4), where 

a±(t) = 2y/TQ(-g(t) T S + cos A:). 

To expand Weber functions we assume below that a± ^> 
1, which implies that 



TQ > 



1 



4(1 -s(f)-*) 2 



for small k modes. This condition can be equivalently 
rewritten to 

l-0(t)»*,0(l)/V7Q. 

To proceed, we need the following identities [31]. For 
|arg(s)| < 37r/4 

V m (s) = e- s2 ^s m [l + (9( 5 - 2 )], 
while for 7r/4 < arg(s) < 57r/4 



Keeping only the leading terms in the expansion of the 
Weber functions in ([TT]) and ([T2]h we find that 



«*(*) 



TrrA e -^Q/8 e - i [«±(«)+-TQlna±(*)]/4 



r(l - ir' Q /4) 



-,(A2) 



U±(t) W e — Q/4 e -/4 e i[ai(«)+-Qlna ± (*)]/4^ (A3) 

Having these expressions we get 

«fc*(*K(*) = [1 " ex P (-^q/ 2 )] exp(-ixW) , (A4) 



exp I 



-7TT, 



£/2)exp (»*(*)), (A5) 



where 



X(t) = [ot{t) -a 2 + (t)]/4 + T' Q Ha-(t)/a + (t)]/4. (A6) 

Next we additionally simplify these results by consider- 
ing k <C 7r/4 and 5 <C 1 limits. These assumptions result 
in Tg « 4/c 2 tq and « 45tq[1 — Substituting 
#(£) = 1 — t/TQ into (|A6|) one gets « 4£5. After some 
algebra these expressions allow for writing the decoher- 
ence factor in momentum space as 



F k (t) « 1 -4 sin 1 



(4«) (e 



-2tttq/c 2 



(A7) 



Large momentum expansion for the modes ex- 
cited after crossing the critical point at g c — — 1: 

We assume here that g(t) < —1. The calculations are 
similar as above with the only difference that slow driv- 
ing through the critical point at g c = — 1 excites momen- 
tum modes near k = tt rather then those close to k = 0. 
The small k modes, excited during crossing the critical 
point at g c = 1, are not affected by driving through the 
critical point at g c = -1: (|Mj) - (1X7)) hold for g(t) < -1. 
The solutions ()A2|) and ()A3|) are valid again, but the 
expansion now requires 



T Q > 



1 



4(-l -g(t)-S) 2 ' 



because we need to substitute cos k ~ — 1 into a±. Equa- 
tions ()A4|) and (jA5|) are the same, but the phase x and 
modified quench time scale t'q have different expansions 
for S < 1 and 7r - fe <C 7r/4: ~ -4tq<J(1 + #) and 

tq ~4(/c-7t) 2 tq. 
Wrapping it up we obtain 



©m(«) = e"^X> m (-s) 



T(-m) 



F fe (i) « 1 - 4 ^-z^oC*-') 2 - e -^r Q (k-Kf j gin 2 ^ 
e *('»+i)T/22>_ TO _ 1 (-i«). where »y(t) = 4r Q £(l + g(t)). 
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